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Abstract 

To make sense of a global space-time model and to give a meaning to 
the coordinates that we use, a choice of a constant curvature space-metric 
of reference it is as much necessary as it is a choice of units of mass, 
length and time. The choice we make leads to contradict the belief that 
the exterior domain of a Static Spherically Symmetric (SSS) space-time 
model of finite radius R depends only on the active mass M a of the source. 
In fact it depends on two parameters M a and a new one Q. We prove 
that both can be calculated as volume integrals extended over the whole 
space. 

We integrate Einstein's equations numerically in two simple cases: as- 
suming either that the source of perfect fluid has constant proper density 
or that the pressure depends linearly on the proper density. We confirm 
a preceding paper showing that very compact objects can have active 
masses M a much greater than their proper masses M p , and we conjecture 
that the mass point Fock's model can be understood as the limit of a 
sequence of compact models when both Q and its radius shrink to zero 
and the pressure equals the density. . 



1 Global SSS space-time models 

Using Weyl's like decomposition, and obvious notations, we consider the line- 
elementQ: 

ds 2 = -A 2 dt 2 + A- 2 ds 2 , A = A(r) (1) 

with: 

ds 2 = B 2 dr 2 + BCr 2 dn 2 , B — B(r), C — C(r) (2) 

Reference space-model. We shall refer to this 3-dimensional metric as the 
space-model, and to restrict it as well as the coordinates used we shall require 
the condition: 

C'=-(B-C) (3) 
Introducing the Euclidean reference metric: 

d§ 2 =dr 2 +r 2 dQ 2 , (4) 
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1 c= 1, G = 1 are used throughout 



we can write the equation above as follows: 

(T% - t) k y k = 0, i,j,k = 1,2,3 (5) 

where f and T are the connection symbols of © and @ in which case this 
equation becomes a restriction on the model but remains true whatever space 
coordinates we use. In particular Cartesian coordinates of (0]) become harmonic 
coordinates of (|2|) and this implies, as it is well known, that they are also 
harmonic coordinates of (JT|). 

Einstein's tensor.- We write below the non identically zero components of 
the Einstein's tensor: 
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where a prime means a derivative with respect to r. They have been somewhat 
simplified eliminating the derivatives of C using ([3]) and its derivative: 

C" = 4(3C-3B + rS') (9) 

The global model- the global perfect fluid models that we consider will be 
solutions of Einstein's equations: 

a: S$ = 8irp, b:Sl = -8ttP, c : S% = £f = -8nP (10) 

A more convenient system of differential equations equivalent to the above 
one is the following: 

2AA" 2A' 2 AAA' 

= 8ir(p + 3P) (11) 
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-8irP (12) 



P' = ~(p + P) (13) 

C' = -(B-C) (14) 
r 

Eq. (jTTJ) is a linear combination of Eqs. a and c. Eq. (fT2"]) is Eq. c left 
unchanged. Eq. (fT3"|) is known to be equivalent to Eq. b as a consequence 
of the conservation equations satisfied by the Einstein's tensor provided that 
S{ = -8nP for r = . Eq. flU) is the same as Eq. © 

As we shall see it is important to realize that Eq. (fTTj) can be written as: 

AU = Att( p + 3P)A- 2 , U = \n(A) (15) 
where A is the Laplace operator of the space- model @: 
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MJ = —d^Vgg^djU), g = det \g id | 



(16) 



or explicitly: 



dU 



B 2 Cr' 2 dr \ dr 



(17) 



2 Boundary conditions 

To define a particular class of SSS models demands to be specific about the 
boundary conditions. This requires i) To state regularity conditions at the origin 
r = 0. ii) To implement Lichnerowicz's continuity conditions at the boundary 
of the source, r = R, that we assume to be compact. And iii) to guarantee 
that from r = R onwards p and P are zero and the solution has the correct 
asymptotic behavior. 

Regularity at the origin - To have continuity at the origin requires to have: 

A' =0, B' o = 0, C" o = 0. (18) 
and taking into account ©there follows that we have to have also: 

Bo = Co (19) 

This means in particular that in a neighborhood of the origin we can write: 



A = Ao + -A 2 r 2 + 0(r 3 ), B = B + -B 2 r 2 + 0(r 3 ), C = B + -C 2 r 2 + 0(r 3 ) 

(20) 

with C 2 = 1/2B 2 to satisfy Eq. ([13]). It follows then from © and (J5J that at 
the origin we always have: 

(Sl)o = (S|)o = (Sf)o = -~^!r (2i) 

that guarantees that the solutions of ([TT]) - (fT4]) will satisfy Eq. b of (fT0|) . 

Matching conditions at r — R - The models will be matched in the sense of 
Lichnerowicz, requiring the continuity of A, B and C and its first derivatives at 
r = R, R being the radius of a compact source. 

Asymptotic conditions: they will be the general asymptotic conditions of 
Schwarzschild's model compatible with ([3]). They were derived in [5], and shall 
be used here 0: 



A-i-^ + i^-i^ + op/O (22) 



B = l-\% + 0(llr>) (23) 
3 r 6 

M 2 4 O 

C=l-^- + ^ + 0(l/r 4 ) (24) 



in this paper is QM^ in [5] 
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where: 



M a = limit r _ ! . 0o AV 2 (25) 
is by definition the active mass of the source, and: 

1 



Q = KmiW 0o -BV 4 (26) 
is a new parameter that has already been considered in other contexts 0-[S] 



3 Integral expressions of M a and Q 

We define the proper mass of the source by: 

/>oo 

M p = 4tt / pr 2 dr (27) 
Jo 

Multiplying both members of (|T5|) by B 2 Cr 2 and integrating from to oo 
we get: 

that using (|2"2")l and (IM1) becomes the well-known Tolman's formula [5]: 

M a = 4tt / (p + 3P)yT 2 BCY 2 dr (29) 
Jo 

On the other hand, integrating by parts the first two terms of <S*Qr 2 over all 
space, and taking into account the conditions stated above we obtain: 

/•OO 

M a = M p + 4tt / a p r 2 dr (30) 
Jo 

where: 

5A' 2 6AA'B' AAA' 7A 2 B' 2 2A 2 B' AAA' A 2 A 2 3A 2 B' 



P B 2 B 3 rB 2 4B 4 rB 3 rBC BCr 2 (J2 r 2 B 2(J r 

( 31 ) 

Similarly integrating by parts the first term of (S 2 — S\)r A and using the 
boundary conditions above we obtain: 



Q = / <r p r z dr (32) 
Jo 

where: 



2r 2 A' 2 2A 2 2r 2 A 2 B' 2 2A 2 r 2 AA'B' 2rA 2 B' 

a P ~ B 2 ~ BC ^4 + Q2 + B 3 + B 3 ( 33 ) 
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4 Numerical models: Examples 



We consider two cases based on the following relationships between the density 
p and the pressure P. 

• Case I: we assume that the density p is constant and positive from r = 
to r — R, where P = 0; 

• Case II: we assume that P = kp — rj from r = to r = R, where P is a 
small fraction of an assumed Pq value at the origin. fc = l/3orfc = l and 
•q is a positive constant included to guarantee that the pressure reaches 
the value zero at some radius R; 

The initial initial conditions will be: 

A' = B' = Co = B =C = 1, P = (8n)- 1 (34) 

compatible with the regularity conditions discussed above. 
The following coordinate transformations: 

t = pi, r = vf (35) 
modifies A, B and C as follows: 

A = pA, B = vfiB, C = vpC, (36) 

but leaves invariant ([3]) as well as ([S])-© and therefore p and P. A covariant 
property that will be used below. 

The integration proceeds in two steps and two runs. First of all we integrate 
Eqs. (ITT1) - (fT4)) until P reaches a chosen estimated small value at some r = R. 
Beyond that we set p — 0, P — and the integration proceeds until the quantity 
e = rAp reaches a chosen estimated small value. Let A € and Let B e be the 
corresponding values of A and B. 

The second round consists in integrating in two steps the same differential 
Eqs. (fTT]) - (fl"4"l) with the same initial values of po and Pq but with the following 
different initial conditions: 

s °=£ (37) 

which is one of the coordinate transformations considered in (1381). with: 



fi = A- 1 , A = A e B^ x (38) 

Examples - The parameter rj in Case II is 0.01 for k = 1/3 and 0.10 for k = 1, 
and Eqs. (ITl"l) - (fl"4"|) have been integrated until P became of the order of 1.0 e — 8. 
The following tables where A is the compactness parameter: 

M 

A = A -f (39) 
summarize some of the informative results: 
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Case II: 
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Appendix 

This appendix is an update of our reference [7] that is meant to prove that 
Fock's model [5] of the Schwarzschild's exterior solution can be considered as 
a legitimate limit model of regular SSS models when the proper volume that 
contains the proper mass shrinks to zero. 
For this model we have: 
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and evaluating (|3~Tj) and f|33[) we obtain: 
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and therefore using (|30l) and 
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we get: 

M„ = -.!/„. Q = 0. 
The first of these formulas suggest that the source of Fock's model is: 

p = hvl a 5ir - M a ) 
where the Dirac density will be denned by the two properties: 

f-OC 

r > Ma S(r - M a ) = 0, and 4tt / 5(r - M a ) 

Jm„ 



dr = l 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



Tolman's formula (|2"9"|) for general regular models and Eq. (|44p above suggest 
that if the proper density for Fock's model contributes to only one third of the 
active mass then the pressure must be: 



P = \M a 5{r - M a ) 



(46) 
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so that: 



P = p and p + 3P = M a S(r ~ M a ) (47) 

Now, let us take a look to the geometry @ corresponding to the coefficients 
B and C above. We have: 



y/det g r = B 2 C sin 9 = (r 2 - Ml) sin 9 (48) 

that becomes zero when r = M a . This means that the system of coordinates 
whose meaning is that derived from the metric of reference (j4|) is not to be 
trusted when discussing properties of the model in a neighborhood of r = M a . 

What we do below is more than we need to discuss this point but we prefer 
here to deal with a global condition instead of a local one as we did in [7] and 
introduce a system of coordinates: 



yi = u(r) sin#cos</>, yi — u(r) sin sin 0, j/3 = u(r) cos 9 (49) 
u being a function of r such that: 



^det g y = 1« y/ det g u = u 2 (50) 

Using the polar form of the condition, a short calculation proves that the func- 
tion u will be given by: 

u = (r 3 - ZMlr + 2M% ) 1/3 (51) 

so that: 



u(M a ) = (52) 

The other way around, to get the inverse function r(u) we need to solve the 
third degree algebraic equation: 

r 3 - ZMl r + 2Ml - u 3 = (53) 

The discriminant of this equation is positive if u > 2 2 / 3 M a and negative other- 
wise down to u = 0. In the first case there exist only one real solution, namely: 



r = ( M 3 + l -u 3 + Iv/-4A/3 U 3 + U 6^ /3 +(~M 3 + l -u 3 - ^-4M^u 3 + u^ " 

(54) 

that has an asymptote r = u, while in the second case there are three real 
functions. But there is only one for which r is a positive increasing function of 
u, namely: 

r = 2M a sin (-n + - arccos ( 1 - --^r^ (55) 
V 6 3 \ 2M* V ' 



whose behavior near u = is: 



1 V3V/ 2 1 



II 



3 



^ M " + 3W^^ +0(u9/2) (56) 
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The two solutions here considered joint smoothly at u = 2 2 / 3 M a . 

Using the radial coordinate u instead of r the space-model © becomes: 

ds 2 = Bldu 2 + B'^cin 2 (57) 

where: 

j,, cos ( \-k + | arccos ( 1 — -■|£f ) ) u 2 

Bu = £ = — — v ; JJ (58) 

and B 2 can be approximated in a neighborhood of u=0 by: 
while U u and dll u /du can be approximated by: 



and: 

d[4_ 3u lV3^, c>((fJj (()]) 



du 4 M a 6 Ma 3/2 

Let us now consider again Eq. ([TT]) that becomes now, for Fock's model, the 
Laplace equation: 

— — ( U ^B- 2 — 
u 2 du \ u du 



W» = ~[« l B?=£.)=0 (62) 



that holds for u > but it is actually undefined for u = 0. 
Let us define U u as the linear functional: 

/oo 
U u tpu du (63) 

where y> is any function such that itself and all its derivatives become zero 
beyond some finite value of u. As a distribution we shall have then: 

/oo 
U u Aipu 2 du (64) 



or: 



< AU u ,tp >=4tt limits ^ U u ^- ^ 2 B,; 2 ^j du (65) 
Using now this particular case of Green's formula: 



TT d ( 2 u -2 d{ f\ d ( 2 u -2 dU u\ d ( 2 2 dip \ d ( 2 2 dU u ^ 

(66) 
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and taking into account (|65[) we have: 

< AU u ,tp >= 47rlimit u ^o (u 2 B- 2 (u u ^f- ~ tp^ I j ((»7) 



du du 



Now since: 



limit^o^ 2 ^ 2 ^) - and limit„^ (^n 2 ^) = M a (68) 
we conclude that: 

< AU U , tp >= 47rM Q ^(0), or AU U = A-k^MJ{u) (69) 
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